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Abstract 

An ubiquitous property of plasmas is the so-called Debye shielding of the electrostatic potential. 
Important aspects of Debye screening concern, in particular, the investigation of non-linear charge 
screening effects taking place in strongly-coupled plasmas, that imply a reduction of the effective 
charge characterizing the Debye-Hiickel potential. These effects are particularly relevant in dusty 
plasmas which are characterized by high-Z particles. The investigation of the effective interactions 
of these particles has attracted interest in recent years especially for numerical simulations. In this 
work we intend to analyze the consistency of the traditional mathematical model for the Debye 
screening. In particular, we intend to prove that the 3D Poisson equation involved in the DH model 
does not admit strong solutions. For this purpose a modified model is proposed which takes into 
account the effect of local plasma sheath (i.e., the local domain near test particles where the plasma 
must be considered discrete). Basic consequences of the model are discussed, which concern the 
asymptotic properties of the solutions determined both for weakly and strongly coupled plasmas. 
As an application the charge screening effect in strongly coupled plasmas is investigated and an 
explicit expression of the effective charge for the asymptotic DH potential is determined. 

PACS numbers: 51.50+v, 52.20-j, 52.27.Gr 
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I. INTRODUCTION 



A basic aspect of plasma physics is the so-called Debye shielding of the electrostatic po- 
tential. This consists in the property of plasmas (or electrolytes |l|), either quasi-neutral 
or non-neutral, to shield the electrostatic field produced by charged particles (to be also 
denoted as test particles) immersed in the same system. This result has fundamental con- 
sequences on plasma phenomenology, since it actually limits the range of static Coulomb 
interactions inside the Debye sphere, i.e., at a distance p < Xd from the test particle, Ad 
being the Debye length. As usual here X D = X D \*j , where the sum is carried out 



on all plasma species and Xd s = y ^z^n ' ^ anc ^ ^ os being respectively the s— species 
temperature and number density (the latter defined in the absence of test particles). In 
fact, when both the particles and the plasma are assumed non-relativistic, small-amplitude, 
stationary (or slowly time- and space- varying) , electrostatic perturbations generated by iso- 
lated test particles, result effectively shielded in the external domain, i.e., at distances larger 
than the Debye length Xd- The renewed interest in this problem [3, S, S , f], Q, [si, [9] is 
particularly related to dusty plasmas or colloidal suspensions which are characterized by 
the presence of a large fraction of highly charged particles (grains), i.e., having an electric 
charge Z d e with \Z d \ ;§> 1. 

In this work we intend to analyze the consistency of the traditional mathematical model 
for the so-called Debye screening problem (DSP) originally formulated by Debye and Hiickel 
(DH model [l|). In particular, we intend to prove that the 3D Poisson equation involved 
in the DH model does not admit physically acceptable solutions, i.e., solutions which are 
provided by ordinary functions and are at least continuous in the domain of existence, i.e., 
are so-called classical (or strong) solutions. For this purpose a modified model is proposed 
which takes into account the effect of local plasma sheath (i.e., the local domain near test 
particles where the plasma must be considered discrete). Basic consequences of the model 
are discussed, which concern the asymptotic properties of the solutions determined both for 
weakly and strongly-coupled plasmas. 

Despite previous attempts to construct approximate or exact solutions to the DH model 



id, 
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121 ]. the related mathematical model appears incomplete and can be shown to be 



physically unacceptable, due to the neglect of the local plasma sheath. In fact, it is obvious 
that sufficiently close to the point-particle the so-called weak-field approximation is violated 
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making the DH model invalid [8|. In the past [13J it was pointed out that in such a case 
the test particle does not produce any electric field, but only complete charge neutralization 
by the plasma, thus producing a Debye length which effectively vanishes. Other objections 
concerned the asserted indeterminacy of the solution for x = due to its divergence in the 
same point 14]. These issues were later addressed in a more general context 15], including 
the 2D case where complete neutralization cannot be achieved. To recover the correct 
physical picture the effect of local plasma sheath must be included. Nevertheless, for suitably 
dense plasmas or in the case of plasma species characterized by very high electric charges 
(high-Z), such as dusty plasmas, the weak-field approximation may be locally violated. 
This circumstance, when the effect of finite local plasma sheath is included, occurs if the 
normalized electrostatic potential results of order unit or larger on the boundary of the 
plasma sheath (produced by at least one of the s plasma species), namely for p = p os . 

It is well-known that, in general, the Debye effect occurs provided suitable physical as- 
sumptions are introduced. In particular, the plasma must be assumed appropriately close to 
kinetic Maxwellian equilibrium, in which each particle species is described by a Maxwellian 
kinetic distribution function carrying finite fluid fields [defined respectively by the number 
density, temperature and flow velocity (N, T, V)] . In the absence of test particles these fluid 
fields must be assumed slowly varying in a suitable sense, or constant, both with respect 
to position (r) and time (i). In this regard it is important to remark that the appropriate 
treatment of the plasma sheath surrounding each test particle is essential also for the va- 
lidity of the mathematical model for the Debye screening problem, i.e., for the existence of 
classical solutions of the Debye screening problem, which do not exist when letting x Q = 



Another significant aspect concerns the issue of the absorption of plasma particles by the 



test partic 



species [16 



e, w 



. ,»hich effectively modifies the local charge density of the background plasma 
,0, JiJ- Since the particle capture mechanism is a manifestly charge- dependent 
and velocity-dependent phenomenon (in particular it depends on the angular momentum of 
the incoming particle), it is obvious that in principle it can produce deviations from local 



Maxwellian equilibrium 



m 



201 ] . However, this phenomenon is expected to become relevant 
only if the radii of the test particle and of the surrounding plasma sheath (p p and p Q ) are 
comparable, i.e., p v j p a ~ 1. Instead, is results negligible when p p / p a <C 1. Since dusty and 
colloidal plasmas are characterized by typical grain size p p smaller than 10~ 4 — 10~ 5 cm and 
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radius of plasma sheath p Q of the order of 10~ 2 — 10~ 3 cm, these effects will be considered 
negligible. 

Goal of this work is the analysis of DSP and the definition of a suitably modified math- 
ematical model to take into account the effect of local plasma sheaths in quasi-neutral 
plasmas. In particular, in Sec. 2 a modified Debye screening problem (modified DSP) is 
presented. In Sec. 3 the basic mathematical result is presented which concerns the non- 
existence of classical solutions of DSP. The proof is reached by noting that DSP can be 
obtained as limit problem obtained from the modified DSP. Basic feature of the approach 
is the representation of the Poisson equation in integral form. This permits to analyze the 
asymptotic properties of the solutions of the modified problem in the limit x os — > + . It is 
found, that the limit solution of the modified DSP for x os —>■ + is a distribution which van- 
ishes identically for all p > and is discontinuous in p = 0. Hence, the limit function <&(x) is 
not an strong solution of the DSP equation. This is therefore a characteristic property of the 
DH model. In particular, as a basic consequence, the effective charge of the DH asymptotic 
solution c vanishes identically in such a limit and results independent of the charge of the 
test particle. 

II. THE MODIFIED DSP 

The traditional formulation of the DSP, based on the Debye-Hiickel model []J regards 
the test particles as point-like and having a spherically-symmetric charge distribution while 
ignoring the effect of local plasma sheath. This implies, from the physical standpoint, 
to neglect the discrete nature of the plasma. Here we shall consider a modified modified 

n 

Debye-Hiickel model, based on the introduction of the notion of local plasma sheath |8|. 
In the sequel we shall consider for simplicity of notation the case of a two species-plasma, 
formed by electrons and Hydrogen ions, having an unique plasma sheath. Thus, we shall 
assume that the test particle is represented by a spherically symmetric charge of radius 
Pp. For a particle in which p p < p the plasma sheath is represented by the spherical shell 
centered at the position (center) of the test particle for which p p < p < p , in which the 
plasma charge density (except for the presence of the test particle) results negligible. In the 
sequel we can also let in particular p p = (point-like test particle) or p p = p Q (finite-size 
test particle). The customary DH model is thus recovered letting p p = and taking the 
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limit p — > (or in dimensionless variables, requiring x p = p p /\d = and x D = p /\d 
— > 0). Denoting $a; o (x) the solution of the Poisson equation, here we intend to determine 
its asymptotic properties in the limi t x p — > + , while also letting x p = (see Lemma). As a 



consequence and in agreement with 



m 



151 ]. in such a case it follows that the limit function 
lim a . o ^ + $ Xo (x) = $(x) vanishes identically for x > 0, i.e., lim a , o ^ + ^oO^) — O.In addition, 
in the same set we intend to prove the identity 

/•X 

(3 - lim / dx'x' 2 sinh $ Xo (x')Q(x' - x Q ) = 0, (1) 

X °^ 0+ Jxo 

where 6( the (weak) Heaviside function. In detail the relevant equations valid 

in each subdomain for the normalized electrostatic potential & Xo ( follows. In 

the internal domain < x < x p the electrostatic potential is assumed constant, namely 
®x (x) =^x (x p ). In the plasma sheath x p < x < x p , § Xo (x) satisfies the customary Poisson 
equation in the presence of the charge density produced by a finite-size spherically-symmetric 
charge 

V x^ = -^(x-x p ). (2) 

Finally, in the external domain x > x a there holds the Poisson equation in the presence of 
the plasma charge density: 

l$ Xo = Q{x-x )sinh$ Xo . (3) 

The boundary conditions, imposed respectively at infinity and at the boundary of the plasma 
sheath, are specified as follows 



X— "OO 



X 2 ^-®X (X 

ax 



lim$, o (x) = 0, (4) 
= -P. (5) 



We notice that, if x p < x Q (for example, x p = 0), $ Xo (^) results by assumption at least 
of class (W{ Xp y) , where K{ Xp } = ]^ P ,oo[ . Here x a , (3 are both assumed constant and 
strictly positive real numbers. The problem defined by (J3]),([2]), together with the boundary 
conditions (jlD,©, will be here denoted as modified DSP. From the physical standpoint 
Eqs.((3]),([2]) may be viewed as the Poisson equation for a spherical ideally conducting charge, 
or for a point particle in the presence of a plasma sheath, of radius r Q (i.e., x Q = r j\D in non- 
dimensional variables) which is in electrostatic equilibrium and is immersed in a spatially 
uniform quasi-neutral and Maxwellian plasma. As for the previous DSP equation, it follows 



that, for solutions satisfying the boundary conditions (JIJ,©, in the domain x G ffi{ Xo } Eq.(j3]) 
can be cast in the integral form 

K(x) = mX ~ Xo) - (6) 

x 

' rx POO 

— / dx'x' 2 + / dx'x' sinh $ x . o (x')G(x' — x Q ). 

. J X J X 

In particular, thanks to continuity at x = x Q of $ x . o (x), one obtains the constraint 

POO 

$ Xo (x ) = T - / dx'x'smh$ Xo {x')Q{x' - x Q ), (7) 

J X 

with r = denoting the Coulomb coupling parameter. It is immediate to establish, the 
existence and uniqueness of Q Xo (x) in the functional class C^°°' ) (]R{ 2:o }), together with its con- 
tinuous dependence on initial data, in particular the continuity with respect to the parameter 
x Q G{o} • Moreover, assuming that the weak- fields approximation applies (this condition is 
manifestly fulfilled identically in the weak-coupling ordering, < T ~ 0(e) 1), and is 
satisfied at least for x 3> 1 suitably large, it is immediate to prove that in this subset an 
asymptotic solution of the modified DSP is provided by the external asymptotic solution 

$, (x) = $if)(x) = -e-^. (8) 

Here denoted as of the modified DSP and c = c(x a , T).is the effective dimensionless charge. 
Hence, $i„ (x) reduces formally to the so-called DH potential when x 3> x Q . In the weak- 
coupling ordering it follows c(x OJ T) = j^r, while for strongly-couple plasmas a lower value 
is expected. Furthermore, it is obvious that the limit function \im Xo ^ + & Xo (x) coincides 
with the solution of DSP, i.e., 

lim $ Xo (x) = $(z). (9) 

x o ^0+ 

III. NON-EXISTENCE OF CLASSICAL SOLUTIONS OF DSP 

Let us now analyze, which consequences can be obtained for the Debye-Hiickel problem, 
formally obtained by letting x Q = in the previous equations [in particular Eq.Q]. This 
requires the knowledge of the asymptotic properties of the solution Q Xo (x) in the limit 
x Q — > + . The following Lemma will be invoked 8]: 
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A. LEMMA - Asymptotic properties of & Xo (x) 

For any strong solution of the modified DSP, <$> Xo (x) obtained letting x p = 0, the limit 
function lim a;o _»o+ &x (x) has the following properties: 

1) There results: 

lim $ Xo {x ) = +00; (10) 

x o ^0+ 

2) for any x > the integral limit (QP is satisfied by Q Xo (x).This implies that the limit 
function = \im Xo ^ + & Xo (x) results such for any x > 0,x G{o} 

$(x) = lim $ Xo (x) = 0. (11) 

x a — »0+ 

3) the following limit is satisfied by the boundary value Q Xo (x ) 

lim x o $ Xo (x o ) = 0. (12) 

x o ^0+ 

4) the limit value of the effective dimensionless charge c(x Q , T) for x Q — > + , obtained 
keeping T finite, is 

lim c(x Q ,r) = 0. (13) 

x o -^0+ 

PROOF 

In fact, as a consequence of the integral equation ([7]) and the continuous dependence of 
Q Xo (x) on the initial data, it follows 

/"OO 

lim x & Xo (x ) = (3 — lim x Q I dx'x' sinh <fr Xo (x')0(x' — x a ), (14) 

x — >0+ x o ^0+ J x 

which implies 

lim $ Xo (x ) = 00, (15) 

x o ^0+ 

poo ^ ^ 

lim / dx'x' sinh $ Xo (x)0(x' — x Q ) = 00, (16) 

*o^0+ J Xo 

i.e., the limit function lim Xo ^o+ &x (x) diverges in x = x a . Therefore, due to the continuity 
of & Xo (x) with respect to a; G [0, oof it follows that infinitesimally close to x } x Q = 0, and 
when infinitesimal of the same order, $ Xo (x) must diverge logarithmically as 

8 - (x) ~ ln {^}- (17) 
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Let us now consider the implications of the integral equation ([6]) for the limit function 

\im Xo ^ + § Xo (x) for arbitrary x G{ Xo } • There follows 

^ 3 1 f x ~ - 

lim $ x (x) = lim / dx'x' 2 sinh $ x . (x')Q(x' - x )— (18) 

x o ^0+ X X x o ^0+ J Xo 

POO ^ ^ 

— / dx'x' lim sinh <& Xo (x')Q(x f — x Q ), 

Jx Xo^0+ 

where, due to the asymptotic estimate (fTT|) . the second term on the r.h.s. necessarily diverges 

1 r x ^ ^ 

lim — / dx'x' 2 sinh $ x „(x')Q(x' — x ) = oo (19) 

c ^o+xJ Xn 



xo~*0+ X 

unless there results for any x ^ x Q , x &{ Xo } 

lim $ Xo ( x ) = 0. (20) 

As a consequence of Eq. (f20~l) . from the integral equation ([6l) it follows necessarily that for all 
x > 0: _ 

lim ^ e ( x Z x °) = lim I rrf x V 2 sinh$ :Bo (x , )e(x , -x ), (21) 

x o ^0+ X x o ^0+ X J Xo 

which proves the limit (1TT) . As a consequence it must result necessarily that the limit 
\im Xo ^ + sinh Q Xo (x) is a Dirac delta. The limit ffT2"|) follows immediately from the boundary 
condition (JTj), while Eq. (1201) implies manifestly the limit ( [131 ). 

As an immediate consequence of ther Lemma it follows that the DSP equation obtained 
letting x a = in Eq.([6]) does not admit classical solutions. 



B. THEOREM - Non-existence of classical solutions of DSP 

In the functional class C^({o}) the DSP problem defined by Eqs.(^)and the boundary 
conditions indicated above [Eqs.(BC-l b)(BC-2 b)] does not admit strong solutions. 
PROOF 

In fact, first, we notice that the limit function 

lim $ Xo (x)=$(x), (22) 

x o ~^0+ 

is manifestly a solution of the DSP equation which satisfies the required boundary conditions 
(jUrefBC-2 b). On the other hand, due to the Lemma, this solution is discontinuous in x = x Q 
and results a distribution. Hence it is not a strong (classical) solution of the modified DSP 
problem. 

8 



The basic implication of the Lemma and the theorem is that the DSP equation, provided 
by the DH model, must be regarded as physically unacceptable, since it does not admit strong 
solutions. In this regard it should be noted that, as a basic principle, physically acceptable 
of solutions of ordinary (or partial) differential equations characterizing the classical theory 
of fields must be suitably smooth strong solutions. The modified Debye screening problem 
here defined, instead, exhibits smooth strong solutions and therefore appears, from this 
viewpoint, consistent. 

IV. CONCLUSIONS 

The essential implication of the present result is that the customary mathematical model 
used for the investigation of the Debye screening problem is incorrect. The correct mathe- 
matical model requires, in fact, the treatment of the local plasma sheath, for which a simple 
model is provided by Eq.([6]). Important physical consequences follow. These concern the 
correct estimate of the Debye screening effect, which occurs close to the local plasma sheath 
and, particularly, for highly-charged test particles immersed in strongly-coupled plasmas. 
In fact, the modified DSP can be used to obtain asymptotic estimates for the effective di- 
mensionless charge c(x ,T) carried by the DH potential in strongly-coupled plasmas js]. 
The resulting charge screening effect appears produced by non-linear effects in the Poisson 
equation. As a consequence, outside the Debye sphere (i.e., for x > 1) the DH potential 
generated by highly charged test particles in strongly-coupled plasmas results strongly re- 
duced with respect to the theoretical value observed in the corresponding weakly-coupled 
systems. 
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Figure captions 



Figure 1 - Comparison between f c (x), c(x ,T) and c^(x 0: T). The data are normalized 
with respect to (3, the normalized charge of the isolated test particle. The figure concerns 
the case with (3 = 1 and x Q = 0.05, yielding T = 20. The horizontal straight line represents 
the asymptotic estimate c < - a ^(x D , T), while the curve below it is the graph of f c {x). It follows 
f c (x ) = 0.564, while the asymptotic value c(x Q , T) = 0.493 is reached approximately at 
x ~ 0.3, and the upper bound for the normalized effective charge is c^(x Q , Y) = 0.589. 

Figure 2 - Comparison between f c (x), c(x Q , T) and c^(x Q , T) for (3 = 5 and x Q = 0.2 (with 
T = 25). In this case f c (x ) = 0.369, while the asymptotic value c(x ,T) = 0.28 is reached 
approximately at x w 0.4, and c (a)( (a; ,r) = 0.38. 

Figure 3 - Comparison between f c (x),c(x ,T) and c (a)( (a; ,r) for (3 = 10 and x Q = 0.3 
(with T = 33). In this case it is found f c (x ) — 0.273, while the asymptotic value c(x ,T) 
= 0.188 is reached approximately at x m 0.5, and c (a)( (x ,r) = 0.303. 
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